home *** CD-ROM | disk | FTP | other *** search
/ Aminet 15 / Aminet 15 - Nov 1996.iso / Aminet / misc / math / Laplace.lha / Laplace / Data / Projects.lha / LinA_Blatt_2_7 next >
Encoding:
Text File  |  1996-07-22  |  823 b   |  60 lines
  1. #LAPLACE
  2. #USEFLOAT=FALSE
  3. #SIMPLIFY=TRUE
  4. #ERRBND=15
  5. #EXPDSP=5
  6. #PREC=8
  7. #F1=lambda
  8. #F2=g(x)=
  9. #F3=h(x)=f(x)/g(x)
  10. #F4=matrix(
  11. #F5=vector(
  12. #F6=diff(f)
  13. #F7=diff(g)
  14. #F8=diff(h)
  15. #F9=taylor(f,0)
  16. #F10=taylor(f,1)
  17. #F11=%(run type s:startup-sequence)
  18. #F12=
  19. #F13=
  20. #F14=
  21. #F15=
  22. #F16=
  23. #F17=
  24. #F18=
  25. #F19=
  26. #F20=
  27. const(lambda)
  28. Lambda=[lambda,0,0;0,lambda,0;0,0,lambda]
  29.  
  30.  
  31. M=[1,-2,2;-2,-2,4;2,4,-2]
  32. Chi(lambda)=-det(M-Lambda)
  33. (lambda-2)*(lambda-2)*(lambda+7)
  34.  
  35.  
  36. M=[1,1/2,1/2;1/2,1,1/2;1/2,1/2,1]
  37. Chi(lambda)=-det(M-Lambda)
  38. (lambda-2)*(lambda-1/2)^2
  39.  
  40. M=[2,2,-2;2,5,-4;-2,-4,5]
  41. Chi(lambda)=-det(M-Lambda)
  42. (lambda-1)^2*(lambda-10)
  43.  
  44. A=[1,-1;-1,2]
  45. B=[-1,2;2,1]
  46. const(a,b,c,d)
  47. P=[a,c;b,d]
  48. trans(P)*A*P
  49. trans(P)*B*P
  50.  
  51. Alpha=eq(a*b-c*b-a*d+2*c*d,0)
  52. Beta=eq(-b*a+2*d*a+2*b*c+d*c,0)
  53.  
  54. b=1,c=1
  55. a=(sqrt(29)-3)/2*c
  56. d=-b*c/(a+3*c)
  57. P=[a,c;b,d]
  58. eval(trans(P)*B*P)
  59.  
  60.